Advanced GNSS-R instruments for altimetric and scatterometric applications

This work is the result of more than eight years during a bachelor thesis, a master thesis, and the Ph.D. thesis dedicated to the development of the Microwave Interferometric Reflectometer (MIR) instrument. It summarizes all the knowledge acquired during this time, and describes the MIR instrument as detailed as possible. MIR is a Global Navigation Satellite System - Reflectometer (GNSS-R), that is, an instrument that uses Global Navigation Satellite System (GNSS) signals scattered on the Earth's surface to retrieve geophysical parameters. These signals are received below the noise level, but since they have been spread in the frequency domain using spread-spectrum techniques, and in particular using the so-called Pseudo Random Noise (PRN) codes, it is still possible to retrieve them because of the large correlation gain achieved. In GNSS-R, two main techniques are used for this purpose: the conventional technique cGNSS-R and the interferometric one iGNSS-R, each with its pros and cons. In the former technique, the reflected signal is cross-correlated against a locally generated clean-replica of the transmitted signal. In the latter technique the reflected signal is cross-correlated with the direct one. Nowadays multiple GNSS systems coexist, transmitting narrow and wide, open and private signals. A comparison between systems, signals, and techniques in fair conditions is necessary. The MIR instrument has been designed as an airborne instrument for that purpose: the instrument has two arrays, an up-looking one, and a down-looking one, each with 19 dual-band antennas in a hexagonal distribution. The instrument is able to form 2 beams at each frequency band (L1/E1, and L5/E5A), which are pointing continuously to the desired satellites taking into account their position, as well as the instrument's position and attitude. The data is sampled and stored for later post-processing. Last but not least, MIR is auto-calibrated using similar signals to the ones transmitted by the GNSS satellites. During the instrument development, the Distance Measurement Equipment/TACtical Air Navigation (DME/TACAN) signals from the Barcelona airport threatened to disrupt the interferometric technique. These signals were also studied, and it was concluded that the use of a mitigation systems were as strongly recommended. The interferometric technique was also affected by the unwanted contribution of other satellites. The impact of these contributions was studied using real data gathered during this Ph.D. thesis. During these 8 years, the instrument was designed, built, tested, and calibrated. A field campaign was carried out in Australia between May 2018 and June 2018 to determine the instrument's accuracy in sensing soil moisture and sea altimetry. This work describes each of these steps in detail and aims to be helpful for those who decide to continue the legacy of this instrument.Este trabajo es el resultado de más de 8 años de doctorado dedicados al desarrollo del instrumento Microwave Interferometric Reflectometer (MIR). Esta tesis resume todo el conocimiento adquirido durante este tiempo, y describe el MIR lo más detalladamente posible. El MIR es un Reflectómetro de señales de Sistemas Globales de Navegación por Satélite (GNSS-R), es decir, es un instrumento que usa señales de GNSS reflejadas en la superficie de la tierra para obtener parámetros geofísicos. Estas señales son recibidas bajo el nivel de ruido, pero dado que han sido ensanchadas en el dominio frecuencial usando técnicas de espectro ensanchado, y en particular usando códigos Pseudo Random Noise (PRN), es todavía posible recibirlas debido a la elevada ganancia de correlación. En GNSS-R existen dos técnicas para este propósito: la convencional (cGNSS-R), y la interferométrica (iGNSS-R), cada una con sus pros y sus contras. En la primera se calcula la correlación cruzada de la señal reflejada y de una réplica generada del código transmitido. En la segunda técnica se calcula la correlación cruzada de la señal reflejada y de la señal directa. Hoy en día muchos sistemas GNSS coexisten, transmitiendo señales de distintos anchos de banda, algunas públicas y otras privadas. Una comparación entre sistemas, señales, y técnicas en condiciones justas es necesaria. El MIR es un instrumento aerotransportado diseñado como para ese propósito: el instrumento tiene dos arrays de antenas, uno apuntando al cielo, y otro apuntando al suelo, cada uno con 19 antenas doble banda en una distribución hexagonal. El instrumento puede formar 2 haces en cada banda frecuencial (L1/E1 y L5/E5A) que apuntan continuamente a los satélites deseados teniendo en cuenta su posición, y la posición y actitud del instrumento. Los datos son guardados para ser procesados posteriormente. Por último pero no menos importante, el MIR se calibra usando señales similares a las transmitidas por los satélites de GNSS. Durante el desarrollo del instrumento, señales del sistema Distance Measuremt Equi Distance Measurement Equipment/TACtical Air Navigation (DME/TACAN) del aeropuerto de Barcelona mostraron ser una amenaza para la técnica interferométrica. Estas señales fueron estudiadas y se concluyó que era encarecidamente recomendado el uso de sistemas de mitigación de interferencias. La técnica interferométrica también se ve afectada por las contribuciones no deseadas de otros satélites, llamado cross-talk. El impacto del cross-talk fue estudiado usando datos reales tomados durante esta tesis doctoral. A lo largo de estos 8 años el instrumento ha sido diseñado, construido, testeado y calibrado. Una campaña de medidas fue llevada a cabo en Australia entre Mayo de 2018 y Junio de 2018 para determinar la capacidad del instrumento para estimar la humedad del terreno y la altura del mar. Este documento describe cada uno de estos pasos al detalle y espera resultar útil para aquellos que decidan continuar con el legado de este instrumento.Postprint (published version

Performance assessment of time–frequency RFI mitigation techniques in microwave radiometry

©2017 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.Radio–frequency interference (RFI) signals are a well-known threat for microwave radiometry (MWR) applications. In order to alleviate this problem, different approaches for RFI detection and mitigation are currently under development. Since RFI signals are man made, they tend to have their power more concentrated in the time–frequency (TF) space as compared to naturally emitted noise. The aim of this paper is to perform an assessment of different TF RFI mitigation techniques in terms of probability of detection, resolution loss (RL), and mitigation performance. In this assessment, six different kinds of RFI signals have been considered: a glitch, a burst of pulses, a wide-band chirp, a narrow-band chirp, a continuous wave, and a wide-band modulation. The results show that the best performance occurs when the transform basis has a similar shape as compared to the RFI signal. For the best case performance, the maximum residual RFI temperature is 14.8 K, and the worst RL is 8.4%. Moreover, the multiresolution Fourier transform technique appears as a good tradeoff solution among all other techniques since it can mitigate all RFI signals under evaluation with a maximum residual RFI temperature of 21 K, and a worst RL of 26.3%. Although the obtained results are still far from an acceptable bias Misplaced < 1 K for MWR applications, there is still work to do in a combined test using the information gathered simultaneously by all mitigation techniques, which could improve the overall performance of RFI mitigation.Peer ReviewedPostprint (author's final draft

Cross-talk statistics and impact in interferometric GNSS-R

©2016 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.This paper presents a statistical analysis of the crosstalk phenomenon in interferometric Global Navigation Satellite Systems Reflectometry (iGNSS-R). Crosstalk occurs when the Delay-Doppler Map (DDM) of a tracked satellite overlaps others fromundesired satellites. This study is performed for ground-based and airborne receivers and for a receiver onboard the International Space Station (ISS) such as the upcoming GNSS Reflectometry, Radio Occultation and Scatterometry onboard ISS experiment. Its impact on ocean altimetry retrievals is analyzed for different antenna arrays. Results show that for elevation angles higher than 60 degrees, crosstalk can be almost permanent from ground, up to 61% from airborne receivers at 2-km height, and up to similar to 10% at the ISS. Crosstalk can only be mitigated using highly directive antennas with narrow beamwidths. Crosstalk impact using a seven-element hexagonal array still induces large errors on ground, but reduces to centimeter level on airborne receivers, and is negligible from the ISS.Peer ReviewedPostprint (author's final draft

SNR degradation in GNSS-R measurements under the effects of radio-frequency interference

©2016 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.Radio-frequency interference (RFI) is a serious threat for systems working with low power signals such as those coming from the global navigation satellite systems (GNSS). The spectral separation coefficient (SSC) is the standard figure of merit to evaluate the signal-to-noise ratio (SNR) degradation due to the RFI. However, an in-depth assessment in the field of GNSS-Reflectometry (GNSS-R) has not been performed yet, and particularly, about which is the influence of the RFI on the so-called delay-Doppler map (DDM). This paper develops a model that evaluates the contribution of intra-/inter-GNSS and external RFI effects to the degradation of the SNR in the DDM for both conventional and interferometric GNSS-R techniques. Moreover, a generalized SSC is defined to account for the effects of nonstationary RFI signals. The results show that highly directive antennas are necessary to avoid interference from other GNSS satellites, whereas mitigation techniques are essential to keep GNSS-R instruments working under external RFI degradation.Peer ReviewedPostprint (author's final draft

Advanced GNSS-R signals processing with GPUs

Global navigation satellite system reflectometry (GNSS-R) is a group of techniques that uses satellite navigation signals as signals of opportunity for remote sensing applications. In GNSS-R, large amounts of data are acquired and have to be processed. Computation time is typically the bottleneck for ground and airborne experiments. This article presents an efficient solution for off-line GNSS-R processing data taking advantage of graphics processing units (GPUs). After comparing to the typically used CPU languages, such as MATLAB and C++, the advantage of using parallel processing on the GPU is clear. GPU-based computation can reduce the processing time by as much as 95% of the acquisition time of the data. An implementation taking advantage of a home-use GPU is proposed for the data processing units. Thanks to its efficiency, even real-time processing experiments are feasible.This work was supported in part bythe Spanish Ministry of Economy and Competitiveness and FEDER EU underthe project “AGORA: Técnicas Avanzadas en Teledetección Aplicada UsandoSeñales GNSS y Otras Señales de Oportunidad” (MINECO/FEDER) ESP2015-70014-C2-1-R by the Agencia Estatal de Investigación; in part by the SpanishMinistry of Science, Innovation, and Universities, “Sensing with Pioneering Opportunistic Techniques,” under Grant RTI2018-099008-B-C21; in part by the Collaboration under Grant 998142 by the Spanish Ministry of Education; inpart by Unidad de Excelencia María de Maeztu MDM-2016-0600; and in partby Prof. Camps ICREA Academia 2015 award of the Generalitat de Catalunya.Peer ReviewedPostprint (published version

Wind direction signatures in GNSS-R observables from space

Peer ReviewedPostprint (published version

A generic level 1 simulator for spaceborne GNSS-R missions and application to GEROS-ISS ocean reflectometry

©2017 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes, creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.In the past decade Global Navigation Satellites System Reflectometry (GNSS-R) has emerged as a new technique for earth remote sensing for various applications, such as ocean altimetry and sea state monitoring. After the success of the GNSS-R demonstrator payloads aboard the UK-DMC or TDS-1 satellites; at present, there are several missions planned to carry GNSS reflectometers. The GNSS rEflectometry, Radio Occultation, and Scatterometry onboard International Space Station (GEROS-ISS) is an innovative ISS experiment exploiting GNSS-R technique to measure key parameters of ocean, land, and ice surfaces. For GEROS-ISS mission, the European Space Agency (ESA) supported the study of GNSS-R assessment of requirements and consolidation of retrieval algorithms (GARCA). For this, it was required to accurately simulate the GEROS-ISS measurements including the whole range of parameters affecting the observation conditions and the instrument, which is called GEROS-SIM. To meet these requirements, the PAU/PARIS end-to-end performance simulator (P$^{2}$EPS) previously developed by UPC BarcelonaTech was used as the baseline building blocks for the level 1 (L1) processor of GEROS-SIM. P$^{2}$EPS is a flexible tool, and is capable of systematically simulating the GNSS-R observations for spaceborne GNSS-R missions. Thanks to the completeness and flexibility, the instrument-to-L1 data module of GEROS-SIM could be implemented by proper modification and update of P$^{2}$EPS. The developed GEROS-SIM was verified and validated in the GARCA study as comparing to the TDS-1 measurements. This paper presents the design, implementation, and results of the GEROS-SIM L1 module in a generic way to be applied to GNSS-R instruments.Peer ReviewedPostprint (author's final draft

Vegetation canopy height retrieval using L1 and L5 airborne GNSS-R

Vegetation canopy height (CH) is one of the important remote-sensing parameters related to forests’ structure, and it can be related to the biomass and the carbon stock. Global navigation satellite system-reflectometry (GNSS-R) has proved capable to retrieve vegetation information at a moderate resolution from space (20–65 km) using L1 C/A signals. In this study, data retrieved by the airborne microwave interferometric reflectometer (MIR) GNSS-R instrument at L1 and L5 are compared to the Global Forest CH product, with a spatial resolution of 30 m. This work analyzes the waveforms (WFs) measured at both bands, and the correlation of the waveform width and the reflectivity values to the CH product. A neural network algorithm is used for the retrieval, showing that the combination of the reflectivity and the waveform width allows to estimate the CH information at a very high resolution, with a root-mean-square error (RMSE) of 4.25 and 4.07 m at L1 and L5, respectively, which is an error about 14% of the actual CH.Postprint (updated version

Parameter considerations for the retrieval of surface soil moisture from spaceborne GNSS-R

The Microwave Interferometric Reflectometer (MIR) is an airborne GNSS-R instrument developed by Universitat Politècnica de Catalunya. In 2018, it was flown twice over the agricultural Yanco area, New South Wales, Australia, once after a very dry period, and a further time the day after a strong rain event. This rain event resulted in many crop fields being entirely flooded, producing a saturation in the GNSS-R reflectivity value. In this work, the received data set is processed to identify the optimum integration time with the goal to minimize pixel blurring. This issue is assessed for airborne conditions, and then extra-polated to the spaceborne case. The presented results show that the blurring of the GNSS waveform is produced even from an airborne sensor with short integration times. Following the determination of an optimal integration time for the platform in use, the surface roughness term in the reflectivity equation can be isolated due to the signal saturation during very wet surface conditions. The final results from the two channels (L1 C/A and L5) are subsequently presented. In this case, it is shown that most reflectivity variations in GNSS-R measurements are linked to surface roughness and Speckle noise fluctuations rather than soil moisture changes.Postprint (updated version

First experimental evidence of wind and swell signatures in L5 GPS and E5A Galileo GNSS-R waveforms

© 2022 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes,creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.As compared to the using L1C/A signals, L5/E5a Global Navigation Satellite System - Reflectometry (GNSS-R), gives improved resolution over the Earth's surface due to the sharper auto-correlation function. Furthermore, the larger transmitted power (+3dB with respect to L1 C/A), and correlation gain (+40dB) allows the reception of weaker reflected signals. If high directivity antennas are used, very short incoherent integration times are needed to have enough signal-to-noise (SNR) ratios, allowing the reception of multiple specular reflection points such as crest of consecutive waves without the blurring induced by long incoherent integration times. This study presents for the first time experimental evidence of the wind and swell waves signatures in the GNSS-R waveforms, and compares them with models.Postprint (author's final draft